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Abstract 

Linear stability of both sign-definite (positive) and sign-indefinite solitary waves near pitch- 
fork bifurcations is analyzed for the generalized nonlinear Schrodinger equations with arbitrary 
forms of nonlinearity and external potentials in arbitrary spatial dimensions. Bifurcations of 
linear-stability eigenvalues associated with pitchfork bifurcations are analytically calculated. It 
is shown that the smooth solution branch switches stability at the bifurcation point. In ad- 
dition, the two bifurcated solution branches and the smooth branch have the opposite (same) 
stability when their power slopes have the same (opposite) sign. One unusual feature on the 
stability of these pitchfork bifurcations is that the smooth and bifurcated solution branches can 
be both stable or both unstable, which contrasts such bifurcations in finite-dimensional dynam- 
ical systems where the smooth and bifurcated branches generally have opposite stability. For 
the special case of positive solitary waves, stronger and more explicit stability results are also 
obtained. It is shown that for positive solitary waves, their linear stability near a bifurcation 
point can be read off directly from their power diagram. Lastly, various numerical examples are 
presented, and the numerical results confirm the analytical predictions both qualitatively and 
quantitatively. 

1 Introduction 

Bifurcation of solitary waves is an important phenomenon in nonlinear wave equations. One type 
of bifurcation is the so-called pitchfork bifurcation, where a smooth branch of solitary waves exists 
on both sides of the bifurcation point, but two additional solution branches bifurcate out to only 
one side of the bifurcation point. The most common pitchfork bifurcation of solitary waves is the 
symmetry-breaking bifurcation, where solitary waves on the smooth branch have certain symmetry, 
but solitary waves on the bifurcated branches lose that symmetry and become asymmetric. This 
symmetry-breaking bifurcation occurs frequently in various nonlinear wave models originating from 
diverse physical disciplines (such as nonlinear optics and Bose-Einstein condensates). For instance, 
this bifurcation has been reported in the nonlinear Schrodinger (NLS) equations with external 
potentials [UEJEISIEIEIEIEIEI [10l [11]. Physically these NLS equations govern nonlinear 
light propagation in refractive-index-modulated optical media \12\ [13] and atomic interaction in 
Bose-Einstein condensates loaded in magnetic or optical traps (in the latter community these 
equations are called the Gross-Pitaevskii equations |14j). This symmetry-breaking bifurcation has 
also been reported in the linearly-coupled NLS equations which govern light transmission in dual- 
core couplers |X5|. I16j. Analytical studies of symmetry-breaking bifurcations have also been made, 
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mostly for the NLS equations with special types of nonlinearities and potentials (see [H HI [6j 
l9[ 111] for instance). In pQ the authors considered a one-dimensional NLS equation with focusing 
cubic nonlinearity and a Dirac-type symmetric double-well potential, and showed the presence of 
symmetry breaking bifurcation as well as the exchange of dynamical stability from the symmetric 
branch to the asymmetric branch at the bifurcation point. In [2] the authors considered a class of 
multi-dimensional NLS equations with focusing cubic nonlinearity and symmetric potentials, and 
showed that symmetry-breaking bifurcation occurs when the power (also called the squared norm in 
mathematics and particle numbers in Bose-Einstein condensation) of the symmetric solitary waves 
increases above a certain threshold, provided that the first two eigenvalues of the linear potential 
are sufficiently close to each other (such as in double-well potentials with large separation between 
the two wells). In addition, the authors showed that above this power threshold, the symmetric 
states become unstable, and a pair of orbitally stable asymmetric states appear. In [6j, the author 
considered a class of multi-dimensional NLS equations with defocusing power nonlinearity and 
symmetric double-well potentials in the semiclassical limit, and showed that symmetry-breaking 
bifurcations occur for antisymmetric solitary waves. In [8], the authors considered a class of one- 
dimensional NLS equations with focusing power nonlinearity and a symmetric potential, and showed 
that symmetry-breaking bifurcation occurs for positive symmetric solitary waves if the potential 
satisfies certain requirements. In addition, they showed that the symmetric branch changes stability 
at the bifurcation point, and the asymmetric branches can be orbitally stable or unstable under 
different conditions. In [9], the authors considered the same class of equations as in [8 J and obtained 
normal forms for these symmetry-breaking bifurcations. In [IT], this author considered the general 
class of NLS equations with arbitrary forms of nonlinearities and potentials in arbitrary spatial 
dimensions, and derived the general analytical conditions for pitchfork bifurcations as well as the 
power formulae for solitary- wave branches near the pitchfork bifurcation point. 

In this paper, we consider the general nonlinear Schrodinger equations with arbitrary forms of 
nonlinearity and potentials in arbitrary spatial dimensions (as in [LT]). These equations include 
the Gross-Pitaevskii equations in Bose-Einstein condensates with attractive or repulsive atomic 
interactions and nonlinear light-transmission equations in linear potentials or nonlinear lattices with 
power or non-power nonlinearities as special cases [HI [131 E] • For this large class of equations, 
we determine the linear stability of both sign-definite (positive) and sign-indefinite solitary waves 
near pitchfork bifurcations. Our strategy is to explicitly calculate the bifurcation of linear-stability 
eigenvalues from the origin, which always takes place whenever a pitchfork bifurcation occurs (see 
Theorem 3 in Sec. 3). Based on this eigenvalue bifurcation and assuming no other instabilities 
interfere, linear stability of solitary waves near pitchfork bifurcations is then obtained (see Theorem 
4 in Sec. 3). We show that the smooth solution branch always switches stability at the bifurcation 
point. In addition, the bifurcated solution branches and the smooth branch have opposite (same) 
stability when their power slopes have the same (opposite) sign. One unusual feature on the linear 
stability of these pitchfork bifurcations is that the smooth and bifurcated solution branches (on 
the same side of the bifurcation point) can be both stable or both unstable, which contrasts such 
bifurcations in finite-dimensional dynamical systems where the smooth and bifurcated branches 
generally have opposite stability [18J. For the special case of positive solitary waves, stronger and 
more explicit stability results are also obtained (see Theorem 5 in Sec. 3). We show that for 
positive solitary waves, their linear stability near a pitchfork bifurcation point can be read off 
directly from their power diagram. Specifically, their linear stability is simply determined by which 
side of the bifurcation point the bifurcated solutions appear and whose power slope of the smooth 
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and bifurcated solutions is larger. Lastly, we present various numerical examples, and show that 
the numerical results confirm the analytical predictions both qualitatively and quantitatively. 

Compared with the earlier analytical results on stability of pitchfork bifurcations (such as in 
[1, 4, 6, 8 ), our stability results have the following three distinctive features. First, our results apply 
to the general NLS equations with no restriction on the nonlinearity, potential or spatial dimensions. 
Second, we made a direct link between linear stability and the power diagram (especially for positive 
solitary waves whose linear stability can be read off entirely from the power diagram). Third, we 
derived explicit analytical formulae for linear-stability eigenvalues of solitary waves, which can be 
useful when quantitative prediction of linear instability is needed. 

2 Preliminaries 

We consider the generalized nonlinear Schrodinger (GNLS) equations with arbitrary forms of non- 
linearity and external potentials in arbitrary spatial dimensions. These equations can be written 
as 

iU t + V 2 U + F(\U\ 2 ,x)U = 0, (2.1) 

where V 2 = d 2 jdx\ + d 2 /dx 2 + • • • + d 2 /dx 2 N is the Laplacian in the iV-dimensional space x = 
(x\,X2, ■ ■ ■ ,xn), and F(-, •) is a general real-valued function which includes nonlinearity as well as 
external potentials. These GNLS equations include the Gross-Pitaevskii equations in Bose-Einstein 
condensates [14J and nonlinear light-transmission equations in linear potentials and nonlinear lat- 
tices [HI [131 E] as special cases. Notice that these equations are conservative and Hamiltonian. 
For a large class of nonlinearities and potentials, this equation admits stationary solitary waves 

U(x,t) = e ilMt u(x), (2.2) 

where u(x) is a real and localized function in the square- integr able functional space which satisfies 
the equation 

V 2 u- fxu + F(u 2 ,x)u = 0, (2.3) 

and fi is a real-valued propagation constant. Examples of such solitary waves can be found in 
numerous books and articles (see [12^13] for instance). In these solitary waves, \i is a free parameter, 
and u(x) depends continuously on fi. Under certain conditions, these solitary waves undergo 
bifurcations at special values of /i. Three major types of bifurcations have been classified |11| . Of 
these bifurcations, stability of solitary waves near saddle-node bifurcations has been analyzed in 
[191 120] . It was shown that no stability switching takes place at a saddle- node bifurcation, which 
dispels a pervasive misconception that such stability switching should occur. In this paper, we 
study the stability of solitary waves near pitchfork bifurcations. 

A pitchfork bifurcation in Eq. ()2. 1 j) is where on one side of the bifurcation point [i = fiQ, there 
is a single solitary wave branch u°(x; /i); but on the other side of ^o, three distinct solitary- wave 
branches appear. One of them is a smooth continuation of the n°(x; /i) branch, but the other two 
branches n ± (x; /i) are new and they bifurcate out at [i = fiQ. 

To present conditions for pitchfork bifurcations, we introduce the linearization operator of Eq. 

Li = V 2 - v + d u [F(u 2 , X )u], (2.4) 
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which is a self-adjoint linear Schrodinger operator. We also introduce the standard inner product 
of functions, 

/oo 
r(x) 5 (x)dx, (2.5) 
-oo 

where the superscript '*' represents complex conjugation. In addition, we define the power of a 
solitary wave it(x; fi) as 

/oo 
u 2 (x;^)dx, (2.6) 
-oo 

and denote the power functions of the smooth and bifurcated solution branches as 

P {fi) = (n (x ;/ u),n (x;^)), P±(/i) = (^(x; //), u±(x; fx)). (2.7) 

If a bifurcation occurs at [i = fJ,o, by denoting the corresponding solitary wave and the L\ 
operator as 

u (x) = n(x;/u ), Lio = il|/Lt=Mo, w=«o> ( 2 - 8 ) 

then Lio should have a discrete zero eigenvalue. This is a necessary condition for all bifurcations, not 
just for pitchfork bifurcations. In the following sufficient conditions for pitchfork bifurcations 
were derived. 

Theorem 1 Assume that zero is a simple discrete eigenvalue of Lio- Denote the real eigen- 
function of this zero eigenvalue as i/j( x )i an d denote 

G(u;x) = F(u 2 ;xK G fc (x) = ^G| u=Uo , = 2,3. (2.9) 

Then if 

(uq,iP) = <G 2 ,^ 3 ) = 0, (2.10) 

fl= (1 - GaL^uo, V 2 ) #0, (2.11) 

and 

5 = <G 3 ,^ 4 ) - 3(G 2 ^ 2 , L$(G 2 ip 2 )) ^ 0, (2.12) 

a pitchfork bifurcation occurs at fi = fio- The new solitary waves ^(x;//) bifurcate to the right 
(left) side of /j, = /j,q if the constants R and S 1 have the same (opposite) sign. 

In addition, slopes of power functions for the smooth and bifurcated solitary-wave branches at 
the pitchfork bifurcation point were also derived in |llj . 

Theorem 2 Suppose the conditions in Theorem 1 hold and a pitchfork bifurcation occurs at 
= no. Then power slopes of the smooth and bifurcated solitary- wave branches at the bifurcation 
point are given as 

P Q (jio) = 2(u ,L^u ), (2.13) 

and 

P> ) = PL(no) = PM + (2.14) 

Here (and in later text) the prime represents the derivative. 

In this article, we consider pitchfork bifurcations in the GNLS equations (|2.ip where the bifur- 
cation conditions in Theorem 1 hold. 
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The main goal of this paper is to determine the linear stability of solitary waves near these 
pitchfork bifurcations. To study this linear stability, we perturb the solitary waves as [13] 

£/(x,i) = e 1 ^ 1 |u(x) + [u(x) + u»(x)]e A * + [u*(x) - u>*(x)]e A *'} , (2.15) 

where v,w <C 1 are normal-mode perturbations, and A is the mode's eigenvalue. Inserting this 
perturbed solution into (12.11) and linearizing, we obtain the following linear eigenvalue problem 



= (2.16) 

where 

£ = f ° L ° 1 * = \ V 1 (2 17) 

L = V 2 -/i + F(u 2 ,x), (2.18) 

and L% is as defined in Eq. (|2.4p . Both Lo and Li are linear Schrddinger operators and are 
Hermitian. In the later text, operator C will be called the linear-stability operator. The eigenvalue 
problem (|2.16p can also be written as 

Lqw = —iXv, L\V = —iXw. (2.19) 

If this linear-stability eigenvalue problem admits eigenvalues A whose real parts are positive, then 
the corresponding normal-mode perturbation in Eq. (I2.15P exponentially grows, hence the solitary 
wave u(x) is linearly unstable. Otherwise it is linearly stable. Notice that eigenvalues of this 
linear-stability problem always appear in quadruples (A, —A, A*, —A*) when A is complex, or in 
pairs (A, —A) when A is real or purely imaginary. 

Using the Lq operator, the solitary wave equation (|2.3p can be written as 



L Q u = 0. (2.20) 

Differentiating this equation with respect to fj,, we find that 

L lUfl = u, (2.21) 

where = du/dji. These two relations will be useful in later analysis. Due to (|2.20p . the linear- 
stability eigenvalue problem (|2.19p admits a zero eigenmode for every solitary wave u(x; fi): 

A = 0, u = 0, w = u. (2.22) 



This zero eigenmode is related to the phase invariance of the GNLS equation (I2.ip , which says that 
if U(x,t) is a solution of ([2.ip . so is U(~x,t)e ie for any real phase constant 9. 

The GNLS equations (|2.ip may be viewed as an infinite-dimensional dynamical system, with 
solitary waves (I2.2p being its fixed points. In this view, it is tempting to deduce the stability 
of pitchfork bifurcations in the GNLS equations (|2.ip from those in finite-dimensional dynamical 
systems. In finite-dimensional dynamical systems, it has been shown that at a pitchfork bifurcation 
point, the smooth fixed-point branch changes its stability. In addition, the two bifurcated fixed- 
point branches have the opposite stability of the smooth fixed-point branch (on the same side of the 
bifurcation point) [18]. However, these stability results were derived under the assumption that zero 
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is a simple eigenvalue of the Jacobian (linearization) matrix of the system at the bifurcation point 
(see Ref. |18j . Theorem 3.4.1, Hypothesis SN1). For the GNLS equations (I2.ip . the counterpart of 
the Jacobian matrix is the linear-stability operator L defined in Eq. (|2,17p . but zero is not a simple 
eigenvalue of C at the bifurcation point (see Eq. (13. 5p below). This means that we cannot apply 
the above stability results from finite-dimensional dynamical systems to pitchfork bifurcations in 
the GNLS equations (12. ip . Instead we have to analyze this stability for Eq. (|2.1j) separately. As 
we will see, stability for pitchfork bifurcations in Eq. (|2,ip shows some novel features which do not 
exist in finite-dimensional dynamical systems. It is relevant to mention that the same situation 
arises for saddle-node bifurcations as well, where it was shown in |19[ [20] that no stability switching 
occurs in the GNLS equations (12. lh even though such stability switching generally takes place in 
finite-dimensional dynamical systems |18j . 



3 Main results 



Our stability analysis starts with the basic fact that, at a pitchfork bifurcation point fi = fiQ, L\q 
has a discrete zero eigenvalue (see earlier text). With the eigenfunction of this zero eigenvalue 
denoted as (see Theorem 1), we have 



Lioip = 0. 



(3.1) 



Thus at the bifurcation point [i = no, in addition to the phase- invariance- induced zero eigenmode 
(|2,22p . the linear-stability eigenvalue problem (|2.19p also admits another bifurcation-induced zero 
eigenmode, 

A = 0, v = ip, w = 0. (3.2) 



Away from the bifurcation point (/i ^ /^o), while the phase-related zero eigenvalue (I2.22P persists 
at the origin, the bifurcation-induced zero eigenvalue (|3.2p bifurcates out since this zero eigenmode 
does not exist in L\ any more. Thus our approach is to determine how this bifurcation-induced 
zero eigenvalue moves out of the origin when [i moves away from fj,Q. We will show that this 
zero eigenvalue only bifurcates out along the real or imaginary axis as a ±A pair. Bifurcation 
along the real axis creates instability, while bifurcation along the imaginary axis does not create 
instability. Thus, based on which direction this zero eigenvalue bifurcates and assuming no other 
instabilities interfere, linear-stability behaviors of solitary waves near the bifurcation point will 
be analytically obtained. In the special case of positive solitary waves, we will show that there 
is indeed no other instabilities interfering near a pitchfork bifurcation, hence stronger and more 
explicit stability results will be derived. 

For later analysis, we introduce two additional notations, 



-^00 — -^o|/i=/^o, «=uo> 



In view of Eq. (|2.20p , we have 

Loo^o = 0, 

thus zero is a discrete eigenvalue of Loo- From Eqs. (12.22P and (I3.2p . we have 



-Co 





Uq 



1> 





0, 



(3.3) 
(3.4) 

(3.5) 
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so zero is also a (multifold) discrete eigenvalue of Cq. 

On the bifurcation of the zero eigenvalue in the linear-stability operator C when \i ^ /io, we 
have the following main results. 

Theorem 3 Assume that zero is a simple discrete eigenvalue of Loo an d Lio- Near a pitchfork 
bifurcation point [i = fj,Q in Theorem 1, if 

(^LooV^O, PM±Q, P±(^o)^0, (3.6) 

then a single pair of non-zero eigenvalues ±A in C bifurcate out along the real or imaginary axis 
from the origin when ^ /iq; 

(a) on the smooth solution branch u°(x; //), the bifurcated eigenvalues A are given asymptotically 
by 

(A ) 2 -» a(n-n ), M^-Mo, (3.7) 

where the real constant a is 

« = T— -prr v ^ o ; (3.8) 

(b) on the two bifurcated solution branches (x; /j,) , the bifurcated eigenvalues are given 
asymptotically by 

(A±) 2 ^/3(/i-^ ), M ^ Mo , (3.9) 

where the real constant j3 is 

P = -2a^- + 0. (3.10) 

Remark 1 Due to the assumption in this theorem, the discrete zero eigenvalue in Cq is not 
embedded inside its continuous spectrum. This fact allows us to calculate eigenvalue bifurcations 
from the origin in C when < \fi — /j,q\ <C 1 by perturbation series expansions (without worrying 
about continuous-wave tails in the eigenfunctions beyond all orders of the perturbation expansion 

nam EaiEg). 

Remark 2 In this theorem, the assumption of zero being a simple discrete eigenvalue of -Lqo 
and Lio is satisfied in all one-dimensional bifurcations and many higher-dimensional bifurcations. 

A direct consequence of Theorem 3 is the following Theorem 4 which summarizes the qualitative 
linear-stability properties of solitary waves near a pitchfork bifurcation point. 

Theorem 4 Suppose at a pitchfork bifurcation point fi = ^o, the solitary wave uo(x) is 
linearly stable (i.e., all its eigenvalues are either zero or purely imaginary); and when \i moves away 
from fj,Q, no complex eigenvalues bifurcate out from non-zero points on the imaginary axis. Then 
under the same conditions of Theorem 3, the smooth solution branch it°(x;//) undergoes stability 
switching at the bifurcation point [with the right (left) side being unstable if the constant a in (|3.8p 
is positive (negative)]. Near the bifurcation point, the two bifurcated solution branches it^x;//) 
and the smooth solution branch (on the same side of the bifurcation point) have opposite (same) 
linear stability when their power slopes Pq(^o) and P±(/io) have the same (opposite) sign. 

Based on this theorem, there are four types of pitchfork bifurcations in the GNLS equations 
(|2.ip . and their schematic diagrams are displayed in Fig. [TJ (note that horizontal flips of these 
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bifurcation diagrams are also admissible, but we will not distinguish them from the ones in Fig. Q] 
for brevity). The bifurcations in Fig. []Ja,b) are qualitatively the same as the supercritical and 
subcritical pitchfork bifurcations in finite-dimensional dynamical systems |18j . In these two cases, 
the bifurcated solution branches and the smooth solution branch have opposite stability. The 
bifurcations in Fig. QJc,d), however, are different. In these two cases, the bifurcated solution 
branches and the smooth branch have the same stability (all stable or all unstable), which seems 
to have no counterpart in the classical dynamical-system theory [18]. Note that the bifurcation in 
Fig. CD|c) has been reported in [8j, but the bifurcation in Fig. QJd) has not been discovered before 
to the author's best knowledge. 




\l \l ll \l 



Figure 1: Schematic diagrams for the four types of pitchfork bifurcations in the GNLS equations 
(|2.ip . Plotted in this figure are deviations u(xo; /u) — uq(xq) versus /i at a representative xo position. 
Solid blue and dashed red lines indicate stable and unstable branches respectively. 

Remark 3 In Theorem 4, the sign of a plays a critical role for the stability outcome. This 
sign can be determined as follows. From formula (|3.8|) . we see that the sign of a is determined by 
the signs of R and (ip,LQQip). We know from Theorem 1 that R and S have the same (opposite) 
sign if the bifurcated solitary waves appear on the right (left) side of /i = ^o- We also know from 
formula (I2.14p that S and P±(no) — Pq(/Uo) have the same sign. Using this information, the sign of 
a can be determined as follows. 

(i) If the bifurcated solitary waves appear on the right side of fi = /xo, then the sign of a is equal 
to the sign of (if) , L^q vp) multiplying the sign of P±(fio) — P^fio); 

(ii) If the bifurcated solitary waves appear on the left side of /u = /Uo, then the sign of a is opposite 
of the sign of (ip , Lqq 1 ip) multiplying the sign of P±({iq) — Pq(hq). 

From this remark, we see that if the sign of (ip, Lqq 1 iP) is known, then the sign of a can be read 
off from the structure of the power diagram (i.e., from which side of [i = fiQ the bifurcated solutions 
reside and which of the power slopes P±(fio) and Pq(//q) is larger). After the sign of a is obtained, 
stability of all the solution branches can then be read off again from the power diagram by Theorem 
4. For instance, when (ip, L^ip) < and the bifurcated solutions appear on the right side of the 
bifurcation point ^ = no, schematic power diagrams of all six possible bifurcation scenarios (with 
stability information indicated) are displayed in Fig. [2j This list of power diagrams is compiled 
according to the signs of Pq(hq), P±(fio), and P4(^o) — -^0(^0)- If (V'i-^oo 1 '^) > 0, these power 
diagrams remain the same, but the stability of all branches is flipped (with "stable" changed to 
"unstable" and vise versa). If the bifurcated solutions appear on the left side of the bifurcation 
point ji = /iQ, their power diagrams (with stability information) can be similarly obtained. 
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Figure 2: Six possible scenarios of the power diagram (with stability indicated) for pitchfork bi- 
furcations in the GNLS equations (|2. 1 j) when (ip, Lqqi[)) < and the bifurcated solutions reside on 
the right side of the bifurcation point. Solid blue and dashed red lines indicate stable and unstable 
branches respectively. 

One might notice that the power diagrams of pitchfork bifurcations in Fig. [2] split out to two 
rather than three branches at the bifurcation point, which is different from the solution-bifurcation 
diagrams in Fig. [TJ The reason is that power slopes of the two bifurcated solution branches u (x; /i) 
are the same at the bifurcation point (see Theorem 2), thus one sees only two power branches instead 
of three [TT] . 

The above results (i.e., Theorems 1-4) are valid for all real- valued solitary waves w(x; /j,) in the 
GNLS equations (|2.ip . including both sign-definite (positive) and sign- indefinite (sign-changing) 
solitary waves. If the solitary waves are positive, then our stability results can be made stronger 
and more explicit. For positive solitary waves in Eq. (|2,ip . it is known that all eigenvalues in the 
linear-stability operator C are either purely real or purely imaginary (see Ref. [13], Theorem 5.2). In 
addition, linear stability of the solitary wave i£o( x ) at the bifurcation point can be determined by the 
generalized Vakhitov-Kolokolov stability criterion [13]. Furthermore, zero is the largest eigenvalue 
of Lq and is simple [21], and (ip, Lqq^) < since operator Loo is semi-negative definite. Using this 
information, together with Theorem 4 and Remark 3, we can obtain the following stronger and 
more explicit theorem for the linear stability of positive solitary waves near a pitchfork bifurcation 
point. This theorem derives linear stability of these solitary waves almost exclusively from their 
power diagram. 

Theorem 5 Suppose solitary waves in the GNLS equations (12. ip are positive near a pitchfork 
bifurcation point fi = fiQ. If Pq(ho) < 0, or zero is the n-th largest discrete eigenvalue of L\q with 
n > 3, then these solitary waves (near the bifurcation point) are linearly unstable. If Pq(hq) > 0, 
P±(fj,o) ^ 0, and zero is the second largest and simple discrete eigenvalue of Lio, then 

(1) when the bifurcated solitary- wave branches u ± (x;//) appear on the right side of \i = no and 
P±([io) < Pq(Mo)) the smooth solution branch u (x;/i) is linearly stable for \i < /j,q and 
unstable for n > fiQ, whereas the bifurcated branches are linearly stable if P^lMo) > and 
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unstable if P±(fio) < 0; 

(2) when the bifurcated solitary waves appear on the right side of fi = /j,q and P4(A*o) > -^0(^0)5 
the smooth solution branch is linearly unstable for fi < fiQ and stable for [i > fiQ, whereas the 
bifurcated branches are always linearly unstable; 

(3) when the bifurcated solitary waves appear on the left side of [i = fio and P±(fio) < -fo(^o)> 
the smooth solution branch is linearly unstable for /i < /io and stable for [i > fiQ, whereas the 
bifurcated branches are linearly stable if P±(fj>o) > and unstable if P±(fio) < 0; 

(4) when the bifurcated solitary waves appear on the left side of /i = fio and P±([io) > Pg(/io), 
the smooth solution branch is linearly stable for /i < jiQ and unstable for [i > fiQ, whereas the 
bifurcated branches are always linearly unstable. 

In terms of solution-bifurcation diagrams, case (1) in this theorem belongs to pitchfork bifur- 
cations of type (a) or (c) in Fig. 1, case (2) belongs to pitchfork bifurcations of type (b) in Fig. 
1, case (3) belongs to pitchfork bifurcations of type (a) or (c) in Fig. 1 with the horizontal axis 
flipped, and case (4) belongs to pitchfork bifurcations of type (b) in Fig. 1 with the horizontal 
axis flipped. Thus for positive solitary waves in the GNLS equations (12, ID . pitchfork bifurcations 
of type (d) in Fig. 1 (and its horizontal flip) cannot occur. 

In terms of power-bifurcation diagrams, case (1) belongs to pitchfork bifurcations of type Fig. 
2(b,c), case (2) belongs to pitchfork bifurcations of type Fig. 2(a), case (3) belongs to pitchfork 
bifurcations of type Fig. 2(d,e) with horizontal-axis flipping as well as branch-stability flipping, 
and case (4) belongs to pitchfork bifurcations of type Fig. 2(f) with horizontal-axis flipping as well 
as branch-stability flipping. 

Remark 4 For positive solitary waves, pitchfork bifurcation cannot occur when the largest 
discrete eigenvalue of L\q is zero. The reason is that for any linear Schrodinger operator, the 
eigenfunction of its largest eigenvalue is always positive (sign-definite) [24]. Thus for this largest 
zero eigenvalue of L10, its eigenfunction ip is positive. Since uo(x) is also positive, then (uo,ip) 7^ 0, 
which violates the conditions of pitchfork bifurcations in Theorem 1 . Thus this case is not mentioned 
in Theorem 5. 

4 Proofs of the main results 

Proof of Theorem 3 The basic idea of the proof is that we first show the algebraic multiplicity 
of the zero eigenvalue in the linear-stability operator £ is four at the bifurcation point [i = fio and 
drops to two away from it, thus a pair of eigenvalues bifurcate out from the origin when fi ^ /xo- 
This pair of eigenvalues must bifurcate along the real or imaginary axis since eigenvalues of C 
would appear as quadruples if this bifurcation were not along these two axes. Then we calculate 
this pair of eigenvalues near the bifurcation point [i = fiQ by perturbation methods. We show that 
the perturbation series for these bifurcated eigenvalues can be constructed to all orders, with the 
leading-order terms given by Eqs. (j3.7j) and (j3.9j) for solitary waves on the smooth and bifurcated 
branches respectively. 

At the pitchfork bifurcation point fi = ^0 hi Theorem 1, (0,uo) T and (^,0) T are eigenfunctions 
of the zero eigenvalue in Cq in view of Eq. (]3.5p . Here the superscript 'T' represents the transpose of 
a vector. Under the assumption in Theorem 3, Cq does not admit any additional eigenfunctions at 
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the zero eigenvalue, thus the geometric multiplicity of this zero eigenvalue is two. Next we determine 
the algebraic multiplicity of this zero eigenvalue by examining its generalized eigenfunctions. 

First, evaluating the relation (|2.2ip along the smooth solution branch n°(x;u) at \i = fj,o, we 

get 

u , (4.1) 



where n° is equal to n° evaluated at fi 



Co 



Thus, 



7/° 







u 



(4.2) 



which means that (u < j l Q,0) T is a generalized eigenfunction of the zero eigenvalue. The next-order 
generalized eigenfunction (f2,g2) T satisfies the equation 



Co 



h 

92 



u 



fj,0 




(4.3) 



so the equation for #2 is 

L o92 = u%. (4.4) 

Since no is a homogeneous solution of this equation, Loo is self-adjoint, and (no, n° ) = Pq(hq)/2 ^ 
by conditions (|3.6p . according to the Fredholm alternative theorem, Eq. (|4.4p does not admit 
localized solutions for g2, thus there are no additional generalized eigenfunctions for (0, no) T - 

Next, we consider generalized eigenfunctions of Co for (^,0) T . The lowest-order generalized 
eigenfunction satisfies the equation 



(4.5) 



According to the assumption in Theorem 3, the kernel of Loo contains a single localized function 
no, and (uQ,ip) = in view of the conditions for pitchfork bifurcations in Theorem 1. Thus from 
the Fredholm alternative theorem, there exists a real localized function L^tp. Consequently, 



" fl ' 






. 91 _ 








Co 





L ooV 







(4.6) 



i.e., {Q,L Q qiI)) t is a generalized eigenfunction of Cq. 
(/2,<?2) T for (V>,0) T satisfies the equation 



The next-order generalized eigenfunction 



Cr 



h 

92 







(4.7) 



so the equation for is 

£10/2 = ^ooV (4-8) 

Since ip is a homogeneous solution of this equation and {ip , Lqq 1 ip) ^ by conditions (13. 6h . Eq. 
(|4.8p does not admit any localized solutions by the Fredholm alternative theorem. Thus there are 
no additional generalized eigenfunctions for (if),0) T . 
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The above analysis shows that Cq has two eigenfunctions and two generalized eigenfunctions at 
the zero eigenvalue, thus the algebraic multiplicity of the zero eigenvalue in Cq (at the bifurcation 
point) is four. 

When // 7^ /j,q, on any solitary-wave branch u°(x;/i) or n ± (x;^), Eqs. (12.20p - (12.2ip hold, thus 
(0, u) T is an eigenfunction of the zero eigenvalue in C (see Eq. (|2.22p ). and (u fl ,0) T is its generalized 
eigenfunction. This eigenfunction and generalized eigenfunction at n ^ fio are the counterparts of 
(0,uo) T and (it° ;0) T at fi = /j,q (see above), and they are induced by the phase invariance of Eq. 
(12. ip . Under conditions in Theorem 3, we can further show, by similar analysis as above, that when 
< |^ — fio\ <^ 1, C does not admit any additional eigenfunctions or generalized eigenfunctions at 
the zero eigenvalue, thus the algebraic multiplicity of this zero eigenvalue in C is two. 

Since the algebraic multiplicity of the zero eigenvalue in C is four at fj, = fio and drops to two 
when fi ^ fiQ, this means that when fi moves away from ^o> a P a h of linear-stability eigenvalues must 
bifurcate out from the origin. Notice that the two multiplicities of the zero eigenvalue associated 
with the phase-invariance mode (|2.22p persist when jj, 7^ hq (see above), it is then clear that the 
two non-zero eigenvalues must bifurcate out from the bifurcation- induced zero eigenmode (13. 2p . 
Since eigenvalues of the linear-stability operator C always appear in quadruples (when they are 
complex) or in pairs (when they are real or purely imaginary) (see discussions below Eq. (I2.19P ). 
this bifurcated pair of eigenvalues then must be real or purely imaginary and be opposite of each 
other as a ±A pair. 

Next, we calculate this pair of bifurcated eigenvalues on the solution branches -u°(x;/i) and 
w ± (x;/i). Since at the bifurcation point the zero eigenvalue of Cq is not embedded inside Cq's 
continuous spectrum, this allows us to calculate this eigenvalue bifurcation by the perturbation 
methods (see Remark 1). 

(a) Eigenvalue bifurcation along the smooth solution branch 

We first calculate this eigenvalue bifurcation along the smooth solution branch u°(x; fi). These 
solitary waves near the bifurcation point n = hq have the following perturbation series expansion 

u°(x;ij) = ^(/i-z^Mx), (4-9) 

k=0 

where U\, u%, . . . are all real functions [llj . As a consequence, operators Lq and L\ on this smooth 
solution branch can be expanded as 

00 00 
L° = Ys(v-M) k L ok , L? = X> - fio) k L lk . (4.10) 

fc=0 k=0 

The linear-stability eigenmodes (v, w, A) bifurcated from the zero eigenmode (13. 2p have the following 
perturbation series expansions: 

00 

W °(x;/,) = ^(/.- M ) fc ^(x), (4.11) 

k=0 

00 

W °(x; u) = A (^ - /i ) 1/2 5> " W))Wx), (4.12) 

AV) = *A (m - mo) 172 ^1 + f> - Mo) fc A^ . (4.13) 
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Below we construct these perturbation series solutions to all orders, and show that the leading-order 
expressions for the bifurcated eigenvalues A°(u) are given by the formula (13. 7p . 

We start by substituting the above perturbation expansions into the linear-stability eigenvalue 
problem (12.19p . From these equations at various orders of \i — fiQ, we get a sequence of linear 
equations for (vk,Wk): 

LiqVq = 0, (4.14) 
L 00 w = v , (4.15) 
L 10 v 1 = X^wq - L u v , (4.16) 
L 00 wi = vi + Xiv - L 01 w , (4.17) 
L 10 v 2 = Xq(wi + Xiw ) - (Lnvi + L 12 v ), (4.18) 
L 00 w 2 = v 2 + Mvi + X 2 v - {L 01 wi + L 02 w ), (4.19) 



71+1 



L\QV n+ i = I Wn + ^2 Afc^n-fc J ~ ^ L lk V n+ i- k , (4.20) 
V fc=l / fc=l 

n+l n+l 

L0QW n+ i = V n+ i + ^kVn+l-k ~ ^2 L 0kW n +l-k, ( 4 -21) 



fc=l k=l 



All these equations are inhomogeneous except the first equation for vq. From the assumption in 
Theorem 3, the v n equations have a single homogeneous solution t/j, and the w n equations have a 
single homogeneous solution no- Since operators Loo and Lio in these equations are self-adjoint, the 
Fredholm alternative theorem says that the inhomogeneous equations above are solvable if and only 
if their right hand sides are orthogonal to their homogeneous solution. These orthogonality condi- 
tions, together with a scaling of the eigenfunction (v,w), will determine the eigenvalue coefficients 
A n as well as functions (v n ,w n ) for all n > 0, as will be demonstrated below. 

First we consider the vo equation ()4.14p . In view of the assumption in Theorem 3, the only 
solution to this equation (after eigenfunction scaling) is 

v = V- (4.22) 

For the wq equation (|4.15p . due to the condition of (no, t/j) = for pitchfork bifurcations in Theorem 
1, the Fredholm condition is satisfied, thus this equation admits a real localized solution Lqq 1 ^, and 
its general solution is 

w Q = LqqV + c n , (4.23) 

where cq is a constant to be determined from the solvability condition of the w\ equation later. 

For the v% equation (I4.16p . it is solvable if and only if its right hand side is orthogonal to ip. 
Utilizing the Vq and wq solutions derived above, this orthogonality condition yields the formula for 
the eigenvalue coefficient Ao as 

A§ = (4.24) 

According to our conditions (|3.6[) . the denominator in the above formula is non-zero, thus Aq is 
well defined and is real. Later in this proof, we will derive a more explicit expression for Aq and 
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show that it is non-zero as well [see Eq. (|4.4ip ]. The above equation (|4.24|) gives two real or purely 
imaginary Ao values as a '±' pair. 

With the eigenvalue coefficient Ao given in (|4.24p , the orthogonality condition of the v\ equation 
(|4.16p is satisfied, thus the general solution for v\ is 

vi = vi + XIcqL^uq + diV>, (4.25) 

where v\ is a real and localized particular solution to the x>\ equation (I4.16P but without the cq term 
in wo on its right hand side (see (|4.23p ). This cq term induces its own particular solution in v\, 
which is the middle term in (|4.25[> . In this term, L^ uq is a real localized function which exists since 
uq is orthogonal to the function tp in the kernel of Lio- The last term in (|4.25p is the homogeneous 
solution, where d\ is a free real constant. Since this homogeneous term in v\ can be lumped to the 
vq term as vq = [1 + di(fi — no)]ip and then eliminated by a scaling of the eigenfunction (v,w), we 
will set d\ = 0. A similar treatment will be applied to all higher v n solutions. 

Now we consider the wi equation (|4.17|) . Its solvability condition is that its right hand side be 
orthogonal to the homogeneous solution uq. Noticing that vq is orthogonal to no and Lq\ is a real 
function, this solvability condition then reduces to 

(uq,V\) - (L m u ,w ) = 0. (4.26) 

To simplify this condition, we recall the relation Lou°(x; u) = 0. By inserting the expansions (14. 9p 
and (|4.10p for Lq and u°(x; u) into this relation and collecting terms of 0(u — /^o), we get 

LqiUq = -LqqUi. (4.27) 

When this relation and the solutions of wq and v\ are substituted into the solvability condition 
(|4,26p and after simple algebra, this solvability condition then reduces to 

c \1(u ,L^qU ) = -(u ,vi) - (u!,ip). (4.28) 

From Eq. f|2. 13|) and conditions (|3.6p . (uq, L^qUq) ^ 0. In addition, we will see from Eq. (|4.4ip 
later that Aq ^ as well. Thus the solvability condition (I4.28P yields a unique real value for cq as 

co = - <y> + <rf . (4.29) 

Consequently, the wq and v\ solutions are now fully determined and are both real. In addition, 
with this Co value, the w\ equation (|4.17p is solvable, and its solution is 

wi = w\ + AiLqqV + ciu , (4.30) 

where w\ is a real and localized particular solution to the w\ equation (|4.17p but without the Ai 
term on its right hand side, and c\ is a constant. The eigenvalue coefficient Ai and the constant c\ 
will be determined from the solvability conditions of the V2 and W2 equations (|4.18p - (|4.19p . 

Next we use the method of induction to show that all higher-order terms in the perturbation 
series expansions (|4.1ip - (|4.13p of (v,w,\) can be successively determined and are all real-valued. 
Suppose the vo,vi, . . . ,v n and u^o, w±, . . . , w n -\ solutions have been fully obtained and are all real. 
In addition, suppose the v n solution is of the form 

V n = V n + Xlc n -iL^Uo, (4.31) 
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where v n is a real and localized function, and c n _i is a real constant. Furthermore, suppose the w n 
solution is of the form 

w n = w n + KLqqiP + c n u , (4.32) 

where w n is a known real and localized function but the coefficients A n and c n are not known yet. 
These assumptions are satisfied when n = 1 (see above), and we now show if they hold for n then 
they would still hold for n + 1 as well. To determine A n , we use the solvability condition of the 
t> n +i equation (14.20p . Inserting (14.32H into this solvability condition, we readily find that 

x Y2t\ LlkVn+l-k) ~ Ag(V>, W n + ELl X kW n - k ) (A qq\ 

which is real. For this A n value, the i> n +i equation (|4.20p is solvable, and its solution is of the form 

v n+ i = v n+ i + \lc n L^u , (4.34) 

where u n +i is a real and localized particular solution to the i> n +i equation but without the c n term 
in w n on its right hand side (see (|4.32|) ). Notice that this v n+ \ solution is of the same form as v n 
in (I4.3ip but with the index n changed to n + 1 . 

To determine the constant c n in the above w n and v n+ \ solutions, we use the solvability condition 
of the w n+ i equation ()4.21 j) . which is that its right hand side be orthogonal to the homogeneous 
solution uq. Inserting the above w n and v n +i solutions into this solvability condition and utilizing 
the relation (|4.2T|) . we find that this solvability condition yields the c n value as 

71+1 

(u , L 01 w n + ^2 L ok w n+1 - k ) 

k=2 



Ao(«o,^i ln o) 

-(u ,v n+1 + ^2\kV n +i-k) - X n (ui,ip) 



k=l 



(4.35) 



which is a real constant. With this c n value, the w n and w n +i solutions are now fully determined 
and are both real. In addition, the w n +\ equation (|4.2ip is now solvable, and its general solution is 

w n+ i = w n+ i + A n+ iLo V + c n+ iu , (4.36) 

where w n +i is a real and localized particular solution to the w n+ i equation ()4.21 j) but without the 
A ra +i term on its right hand side, and c n +i is a constant. This w n +i solution is of the same form 
as w n in (|4.32|) but with the index n changed to n + 1. This completes the induction process. 

It is noted that in the above construction of the bifurcated eigenmodes, while Ao has two 
solutions Ao± (with Ao+ = — Ao-) in view of Eq. (I4.24p . the higher coefficients Ai, A2, • • • as well as 
all (v n ,w n ) functions in the expansions (I4.11|) - (l4.13p depend on Aq only and are thus the same for 
both values of Ao±. It is then clear that the above construction yields two eigenmodes (v±,w±, A^_) 
which correspond to the two choices of the Ao values, and these two eigenmodes are related as 

X° + = -X°_, v° + = v°_, w° + = -w°_. (4.37) 

In addition, v± is always real, and w±, A^_ are either real or purely imaginary. The asymptotic 
formula for the eigenvalues A^_ is 

(A°) 2 ^-Ao(M-/^o), M^Mo, (4.38) 
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where Aq is given in (|4,24p . 

Finally, we simplify the Aq formula (|4.24|) and show that — Aq is equal to a as given in Eq. 
To do so, we expand operator L\ in (|2.4p around fi = fio. Recalling the notations (|2.9p . we readily 
find L\\ in the expansion (14.10P of L\ as 

L n = G 2Ul - 1, (4.39) 

where U\ is the 0(fi — no) term in ti°(x;/i)'s expansion (|4.9p . The expression for u\ has been 
obtained in Ref. [TIj as 

ui = L^Uq + b x ip, (4.40) 

where b\ is a real constant. Inserting the above two expressions into ()4.24p and recalling the 
condition of (G2,^ 3 ) = for pitchfork bifurcations in Theorem 1 as well as the definition of 
constant R in Eq. (|2. 11 [) , we get a more explicit formula for Aq as 

A ° = -TTT^TT- (441) 

In view of the conditions for pitchfork bifurcations in Theorem 1, we see that \l ^ as was 
mentioned before. Inserting this Aq formula into (I4.38p . the final asymptotic expression (13.7P for 
the bifurcated eigenvalues A°(/i) is then derived, with the constant a given by Eq. (|3.8p . 
(b) Eigenvalue bifurcation along the bifurcated solution branches 
Now we calculate eigenvalue bifurcations along the bifurcated solution branches (x; /j,) . For 
convenience, we assume that these bifurcated solutions appear on the right side of [i = (when 
RS > 0, see Theorem 1). The other case of these bifurcated solutions appearing on the left side of 
/i = /io can be similarly treated with trivial modifications, and both cases yield the same eigenvalue 
formula given in Eq. (13. 9|) . 

The bifurcated solitary waves near the bifurcation point [i = [Iq have the following perturbation 
series expansion [11] 

oo 

n ± (x; M ) = J> - W))" 7 V(x). (4.42) 

k=0 

Since these solutions are assumed to exist on the right side of fj, = /j,q, all functions u±, U2, ■ ■ ■ in 
this expansion are real-valued. Operators Lq and L\ on these bifurcated solution branches are 
expanded as 

oo oo 

L^ = Y,^-^) k/2 L k, L± = Y,{n-M) k/2 Liki (4.43) 

fc=0 fc=0 

and all terms in these expansions are real-valued too. Note that quantities u^, Lq^, L^, k = 1, 2, . . . 
in these expansions are different from those in previous expansions (I4,9p - (|4.10p . The linear-stability 
eigenmodes {v , w, A) bifurcated from the zero eigenmode (|3,2p now have the perturbation series 
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expansions 

oo 

^(x; M ) = ^( /U -^ ) fc /V(x), (4.44) 



fc=0 

oo 



x; p) = A - ^ ) fe/2 ^(x), (4.45) 

fc=0 

A ± (//) = »Ao(m - W)) 1/2 U + f> " Mo) fe/ %l . (4.46) 

Below we construct these perturbation series solutions to all orders. 

Before this construction, we first derive a few relations on functions uo,u\,U2 and 113 in (I4.42p . 

which will be needed in later calculations. By inserting expansions ()4.42p and ()4.43p into equations 
(|2.20p and (I2.21j) and at suitable orders, we get the following relations 

£01 «o = -Looui, (4.47) 

^02^0 = -L 01 ui - L ou 2 , (4.48) 

L u ui = 2(u -L10U2), (4.49) 

L\ 2 u\ = 2u\ - 2L\\U2 - 3Li n 3 . (4.50) 

In addition, 

ui = 61V, (4.51) 



where b\ = ±y/6R/S which is non-zero [TT]. Notice that b\ is also real- valued here since RS > 
by our earlier assumption. 

We now substitute the perturbation expansions (|4.44p - (|4.46j) into the linear-stability eigenvalue 

problem (|2.19p . From various orders of (p— /io) 1 / 2 , we get a sequence of linear equations for (v k ,Wk) 
as 

L 10 v = 0, (4.52) 

Loowo = 0, (4.53) 

L 10 vi = \lw - L u vo, (4.54) 

L 00 wi = v - Lqiwo, (4.55) 

L 10 v 2 = \q(wi + Xiw ) - {Lnvi + L 12 vq), (4.56) 

^00^2 = vi + XiVq - {L01W1 + L 2W ), (4.57) 



(n \ n+l 

W n + "^2 ^kW n -k J - ^2 Li k V n+1 - k , (4.58) 
k=l J k=l 

n n+l 

LoO^n+1 = V n + ^ ^kV n -k ~ ^ L ok W n+1 _ k , (4.59) 
k=l k=l 



In view of the assumption in Theorem 3, the solution to the vq equation (|4.52p . after eigenfunction 
rescaling, can be taken as 

v = ui, (4.60) 
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where u% is given in (|4.5ip . The solution to the wq equation (|4.53[) is 

w = c u , (4.61) 

where cq is a constant to be determined. When these (t>o, wq) solutions are inserted into the (v\, w{) 
equations (|4.54j) - (|4.55j) and relations (|4.47|) . (|4.49j) utilized, we find that the solvability conditions 
of the (v±, wi) equations are automatically satisfied due to the condition of {uo,ip) = for pitchfork 
bifurcations in Theorem 1, and these (vi,wi) equations admit localized solutions of the form 

vi = 2(u 2 - L^qU ) + coXqL^qUo, (4.62) 

and 

wi = LqqUi + c ui + ciu , (4.63) 

where L^qUq and LqqU\ are localized functions and c\ is another constant. It is noted that the 
homogeneous solution (in proportion to u\) to the v\ equation is not included in the above v% 
solution since this term can be lumped into the vq term and then eliminated by a rescaling of the 
eigenfunction (v, w) — the same treatment we have applied previously in case (a) (see Eq. (|4.25p ). 

Now we consider the (v 2 ,w 2 ) equations (I4.56|) - (l4.57p . Inserting the above (vo, v\, Wq, wi) so- 
lutions into the right hand side of the v 2 equation, utilizing the relations (j4.47|) - ()4.50|) and after 
simple algebra, the solvability condition of this v 2 equation, which requires that its right hand side 
be orthogonal to the homogeneous solution u\, yields 

^2 /« , 2N (m,ni) +2(u ,u 2 ) ~ 2{u ,L^u ) 

A = (2-c A ) — j — . (4.64) 

(m,L 00 ui) 

It is noted that u\ is proportional to tp and is nonzero, see Eq. (I4.5ip . Thus due to the conditions 
(|3.6p . the denominator in the above equation is nonzero, i.e., (u\, Lq Q ui) 7^ 0. From the expansion 
(I4.42p of the solutions n ± (x;/i) and the condition (uo,ip) = 0, we see that the expansion for the 
power function P±(n) is 

P±(M) = {uo,u ) + P^ )^- f i ) + 0[^-fi f/% (4.65) 

where 

P^o) = (u 1 ,u 1 )+2(u ,u 2 ) (4.66) 

is the power slope at the bifurcation point = fj,Q. From Eq. (I2.14|) in Theorem 2, we also know 
that ^ 

P|( M o) = PL(jm>) = 2(u ,L^ Uo ) + (4.67) 

Using these relations as well as Eq. (|4.51|) , the solvability condition (|4.64p of the v 2 equation reduces 
to 

Ag = (2-c Ag) (4.68) 
iV, L 00 iff) 

From the conditions for pitchfork bifurcations in Theorem 1, R ^ 0, thus Aq 7^ 0. 

Carrying similar calculations to the w 2 equation ()4.57p and recalling the formula (j2. 13|) . the 
solvability condition of this w 2 equation yields 
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When this equation is inserted into (|4.68p . an expression for Aq is then obtained as 

which is real and nonzero. Inserting this Aq formula into (|4.69p . the constant cq can be obtained 
and is also real. 

When Ao and c$ are given by Eqs. (|4.69p - (|4.70p . the solvability conditions of the (^2,^2) 
equations (I4.56|) - (l4.57p are satisfied. Utilizing the relation fj4.47f) . the {v 2 ,w 2 ) solutions are of the 
form 

v 2 = v 2 + Aq(ci + Xico)L^qU , (4.71) 

and 

w 2 = w 2 + Ai-Lqq u\ + c\u\ + c 2 u , (4.72) 

where v 2 and w 2 are localized functions which satisfy equations (|4.56|) - ()4.57[) but without the c\ 
and Ai terms on their right hand sides, and c 2 is another constant. These constants Ai,ci and c 2 
will be determined from the solvability conditions of the higher (v n ,w n ) equations. 

Using the method of induction and after straightforward algebra, we can show that the pertur- 
bation series solution (|4.44|) - ()4.46[) for the eigenmode (v, w, A) can be determined to all orders. In 
addition, the (y n ,w n ) terms for n > 2 are of the form 

v n = v n + Ao(c n _i + A n __ic )Lj" 1 n , n > 2, (4.73) 

and 

w n = w n + Xu-iLqqUx + c n ^iui + c n u Q , n > 2, (4.74) 

where v n and w n are certain localized functions, and A n _i,c n _i,c n are constants. Since this in- 
duction calculation is analogous to that for case (a) in the earlier text, the details are omitted 
here. 

In the above construction of the bifurcated eigenmodes, since Ao has two solutions from Eq. 
(|4.70p . two eigenmodes are then obtained. The eigenvalues of these two modes are either real 
or purely imaginary and are opposite of each other. From the perturbation expansion of these 
eigenvalues in Eq. (|4,46p as well as Eq. (|4.70|) . we see that the asymptotic formula for these 
eigenvalues near the bifurcation point is 

± 2 ^ _2P^o) __R ( } (475) 



which is the same as the formula (|3.9|) in Theorem 3. This completes the proof of Theorem 3. □ 

Proof of Theorem 4 From the assumptions in Theorem 4, the solitary wave uq(x) at the 
bifurcation point fj, = /j,q is linearly stable; and when < \fi — /Xo | "C 1, the only instability- 
inducing eigenvalue bifurcation is from the origin. In Theorem 3, we have shown that from the 
origin, a single pair of eigenvalues ±A in C bifurcate out along the real or imaginary axis. Thus 
the linear stability of these solitary waves near \x = fiQ is determined entirely by whether this pair 
of eigenvalues are real or purely imaginary. On the smooth solitary wave branch ii°(x; fi), this pair 
of eigenvalues are given asymptotically by Eq. (13.7(1 . Thus if a > 0, these eigenvalues are real 
when n > fj,Q and imaginary when fi < /xq, hence the solitary waves are linearly unstable when 
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H > Ho and linearly stable when fj, < Ho- If a < 0, the situation is just the opposite. In both 
cases, stability switches at the bifurcation point. On the bifurcated solitary branches tr^x; /i), the 
bifurcated eigenvalues are given asymptotically by Eq. (|3.9|) . The formula (|3.10|) for the constant 
(3 shows that when the two power slopes Pq(ho) an d P±(W)) have the same sign, (3 and a would 
have the opposite sign, meaning that eigenvalues for the tr^x;//) and u°(x;/i) branches bifurcate 
along perpendicular directions from the origin, hence these solution branches have opposite linear 
stability. On the other hand, if the two power slopes Pq(mo) an d P±(Mo) have the opposite sign, /3 
and a would have the same sign, hence the ^(xj/x) and u°(x;//) branches would have the same 
linear stability. This completes the proof of Theorem 4. □ 

Proof of Theorem 5 If zero is the n-th largest discrete eigenvalue of L\q with n > 3, then L\q 
has two positive eigenvalues, hence the positive solitary wave «o( x ) a t the bifurcation point h = Ho 
is linearly unstable by the generalized Vakhitov-Kolokolov stability criterion |13j . If zero is the 
second-largest discrete eigenvalue of Lio, then L%q has one positive eigenvalue. Meanwhile, by the 
conditions of pitchfork bifurcations in Theorem 1, {uo,ip) = 0. In this case, if Pq(ho) < (> 0), then 
the positive solitary wave ito(x) is linearly unstable (stable) by the generalized Vakhitov-Kolokolov 
stability criterion [13]. When uo(x) is linearly unstable, solitary waves near the bifurcation point 
H = Ho ar e clearly also linearly unstable. 

Next we consider the case when Pq(ho) > 0, Pj_(/io) 7^ and zero is the second largest and 
simple discrete eigenvalue of Liq. In this case, uq(x) at the bifurcation point is linearly stable (see 
above). In addition, for positive solitary waves, all eigenvalues in C are real or purely imaginary 
(see Ref. [13], Theorem 5.2), thus no complex eigenvalues can bifurcate out when \i ^ ho- So the 
assumptions in Theorem 4 are satisfied. Since «o(x) is positive, zero is then the largest eigenvalue 
of Loo an d is simple [23], and (ip, Lqq 1 ^} < as Loo is semi-negative definite and ip ^ 0. In addition, 
zero is a simple eigenvalue of Lio, Pq(ho) an d P±(ho) 7^ by our assumptions above. Thus 
the conditions of Theorem 4 are also met. Hence Theorem 4 can be applied. Using this theorem, 
together with Remark 3 and the fact of {iP,Lqq-iP) < 0, we can then prove the results in the four 
cases of Theorem 5 as below. 

(1) When the bifurcated solitary waves u ± (x; fj,) appear on the right side of fx = pLo and P±(/io) < 
Pq(Ho), a is positive by Remark 3. Then by Theorem 4, the smooth solution branch u°(x; //) 
is stable for fi < Ho an d unstable for fi > Ho- Regarding the bifurcated branches, they and 
the unstable smooth branch (on the right side of h = Mo) should have the opposite (same) 
stability when Pq(ho) an d P±{ho) have the same (opposite) sign. Since Pq(ho) > 0, these 
bifurcated branches are then stable when P±(ho) > an d unstable when P^_(/io) < 0. 

(2) When the bifurcated solitary waves appear on the right side of h = Ho an d P±(ho) > Po (/•*()) > 
0, a is negative by Remark 3. Thus the smooth solution branch is unstable for h < Ho an d 
stable for h > Ho- Since P±(ho) an d Po(ho) are now both positive, the bifurcated branches 
then have the opposite stability of the stable smooth branch (on the right side of h = Ho) are 
thus always unstable. 

(3) When the bifurcated solitary waves appear on the left side of h = Ho an d P'±{ho) < Pq(Mo)) 
a is negative by Remark 3. Then by Theorem 4, the smooth solution branch is unstable for 
H < Ho an d stable for h > Ho- The bifurcated branches (with h < Ho) ar e stable (opposite 
of the unstable smooth branch) if P±(ho) > and unstable (same as the unstable smooth 
branch) if P±{ho) < 0. 



20 



(4) When the bifurcated solitary waves appear on the left side of \i = and P±(no) > Pq(p-o) > 0, 
a is positive by Remark 3; hence by Theorem 4, the smooth solution branch is stable for 
fj, < hq and unstable for fj, > fiQ. Since both P^CA'o) an d Po(fJ-o) are now positive, the 
bifurcated branches are always unstable (opposite of the stable smooth branch on the left 
side of fj, = fj,o). This completes the proof of Theorem 5. □ 

5 Numerical examples 

In this section, we use a few numerical examples to illustrate and confirm the above analytical 
stability results. These examples contain both positive and sign-indefinite solitary waves under 
various self-focusing and self-defocusing nonlinearities in one and higher spatial dimensions. 

Example 1 Our first example is the one-dimensional GNLS equation (|2.ip with a symmetric 
double-well potential and cubic-quintic nonlinearity, 

iU t + U xx - V{x)U + \U\ 2 U - 0.25|C/| 4 C7 = 0, (5.1) 

where the symmetric double-well potential V{x) is taken as 

V(x) = -2.8 [sech 2 (x + 1.5) + sech 2 (x - 1.5)] (5.2) 

and is shown in Fig. 3(a). Pitchfork bifurcations of positive solitary waves in this equation have 
been reported in [IT], and the power diagram for these bifurcations is displayed in Fig. 3(b). This 
power diagram shows that two pitchfork bifurcations occur at points 'A,B' of the power diagram. 
At positions 'c,d,e', profiles of solitary waves are shown in Fig. 3(c,d,e) respectively. It is seen that 
solitary waves at positions 'c,d' are symmetric, while solitary waves at position 'e' are asymmetric. 
Thus symmetry-breaking bifurcations occur at both 'A,B' points. In particular, the c-d power 
branch is the symmetric (smooth) branch, and the e-branch is the asymmetric (bifurcated) branch. 

Since solitary waves in this example are positive, Theorem 5 applies. It is known that for both 
pitchfork bifurcations at points 'A,B', zero is the second largest eigenvalue of L\q [IT] . In addition, 
the power diagram in Fig. 3(b) shows that at both 'A,B' points, < P±(fio) < Pq(^o)- Thus 
Theorem 5 predicts that near the bifurcation point 'A', the smooth (symmetric) branch is linearly 
stable on the left side of 'A' and unstable on the right side of 'A', and the bifurcated (asymmetric) 
branch (which is on the right side of 'A') is linearly stable. Regarding the bifurcation point 'B', 
Theorem 5 predicts that the symmetric branch is unstable on the left side of 'B' and stable on the 
right side of 'B', and the bifurcated asymmetric branch (which is on the left side of 'B') is linearly 
stable. 

Numerically we have found that these analytical predictions are entirely correct. Specifically, 
we have found numerically that the segment of the symmetric branch between points 'A,B' is 
linearly unstable, and the other segments/branches of the power diagram are all linearly stable. 
This stability information is indicated by the solid blue and dashed red lines in Fig. 3(b) for stable 
and unstable parts respectively. To support this stability result, the linear-stability spectra for 
solitary waves at locations 'c,d,e' of the power diagram are numerically computed by the Fourier 
collocation method [13], and the results are displayed in Fig. 3(f,g,h). It is seen that the spectrum 
at 'd' contains a positive eigenvalue, hence its symmetric solitary wave is linearly unstable. The 
spectra at 'c,e', on the other hand, lie entirely on the imaginary axis, thus those solitary waves 
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Figure 3: Pitchfork bifurcations of solitary waves and their linear-stability behaviors in Example 
1. (a) The potential (|5.2|> : (b) the power diagram (solid blue and dashed red indicate stable and 
unstable solutions respectively); (c,d,e) profiles of solitary waves at positions of the same letters in 
(b); the two asymmetric solitary waves in (e) are mirror images of each other with respect to x; 
(f,g,h) linear-stability spectra of the solitary waves in (c,d,e) respectively. 

are linearly stable. These numerical stability results agree completely with the above analytical 
predictions. 

The reader may notice that the asymmetric 'e'-branch in Fig. 3(b) contains two additional 
saddle-node (fold) bifurcations. As was explained in [191 ED], there is no stability switching at 
saddle-node bifurcations in the GNLS equations (12. ip . thus it is not surprising that the entire 
asymmetric 'e'-branch in Fig. 3(b) is linearly stable despite these saddle-node bifurcations. 

Example 2 Our second example is the one-dimensional GNLS equation (|2.1|) with self-focusing 
cubic nonlinearity and a periodic potential, 

iU t + U xx -V(x)U+\U\ 2 U = 0, (5.3) 

where the potential V(x) is 

V(x) = 6sm 2 x. (5.4) 
This equation admits a family of sign-indefinite solitary waves of the form (|2.2p in the semi-infinite 
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bandgap [10J, and its power diagram is shown in Fig. 4(a). This power diagram contains three 
branches which are connected with each other. At the intersection point 'A' between the lower and 
middle branches, a pitchfork bifurcation occurs. This pitchfork bifurcation is better seen in Fig. 
4(b), which shows an amplification of the power diagram in Fig. 4(a) around this intersection point. 
Solitary waves on the lower power branch are anti-symmetric (see Fig. 4(c)), whereas solitary waves 
on the middle power branch are asymmetric (see Fig. 4(d)), hence a symmetry-breaking bifurcation 
occurs at point 'A'. 

At this bifurcation point 'A', we have checked numerically that (ip , Lqq 1 ip) < 0. In addition, the 
power diagram in Fig. 4(b) shows that the bifurcated (asymmetric) solitary waves appear on the 
right side 'A', and .P^Q^o) > Poi^o) > 0- Thus, Remark 3 gives a < 0. We have also checked that 
the solitary wave at point 'A' is linearly stable, and near this point no complex eigenvalues appear. 
Then Theorem 4 predicts that the lower branch is unstable on the left side of 'A' and stable on the 
right side of 'A'. In addition, the bifurcated asymmetric branch is unstable. 
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Figure 4: Pitchfork bifurcations of solitary waves and their linear-stability behaviors in Example 
2. (a) The power diagram (solid blue for stable branches and dashed red for unstable ones); (b) 
amplification of the power diagram in (a) near the pitchfork bifurcation point 'A'; (c,d) profiles 
of solitary waves at locations 'c,d' of the power diagram in (a); (e,f) linear-stability spectra for 
solitary waves in (c,d) respectively. 



Numerically we have found that these analytical predictions are again all correct. Specifically, 
we have determined the stability of these solitary waves through computation of their linear-stability 
spectra, and the stability results are indicated on the power diagram in Fig. 4(a,b), where the stable 
and unstable branches are marked as solid blue and dashed red lines respectively. We see that these 
stability results are in full agreement with the above analytical predictions. To corroborate these 
stability results, the full linear-stability spectra for solitary waves at locations 'c,d' of Fig. 4(a) are 
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displayed in Fig. 4(e,f). These spectra confirm that the anti-symmetric solitary wave [Fig. 4(c)] 
at location 'c' of the lower power branch is indeed linearly stable, whereas the asymmetric solitary 
wave [Fig. 4(d)] at location 'd' of the middle branch is indeed linearly unstable. The unstable 
eigenvalue at location 'd' is real positive, and it bifurcates out from the pitchfork bifurcation point 
'A', in agreement with our theory. 

Example 3 Our third example is the seventh-power GNLS equation with a symmetric double- 
well potential, 

tU t + U xx -V(x)U + \U\ e U = 0, (5.5) 

where the potential V{x) is 

V(x) = -3 sech 2 (x + 1) - 3 sech 2 (x - 1), (5.6) 

which is shown in Fig. 5(b). This equation admits a family of positive solitary waves whose power 
diagram is displayed in Fig. 5(a). This power diagram shows that a pitchfork bifurcation occurs at 
the point 'A'. On the upper b-ci branch, solitary waves are symmetric (see Fig. 5(b,c)), whereas 
on the lower C2 branch, solitary waves are asymmetric (see Fig. 5(c)). 

At the pitchfork bifurcation point 'A', we have checked that zero is the second largest eigenvalue 
of Lio- Thus from the power diagram in Fig. 5(a), Theorem 5 predicts that the symmetric b-ci 
branch is stable on the left side of 'A' and unstable on the right side of 'A'. In addition, the 
asymmetric C2 branch is unstable. These analytical predictions fully agree with the numerical 
stability results in Fig. 5(a) (where the stable and unstable solutions are indicated). These stability 
results are further corroborated in Fig. 5(d,e,f), where linear-stability spectra for solitary waves at 
locations 'b, c\, c% are displayed. 

In this example, after the pitchfork bifurcation occurs (i.e., on the right side of 'A'), both 
the symmetric and asymmetric solution branches are linearly unstable. A similar bifurcation was 
reported numerically in (8] for the eleventh-power nonlinearity but was not found for the present 
seventh-power nonlinearity. 

Example 4 Our last example is the two-dimensional GNLS equation with self-defocusing cubic 
nonlinearity and a symmetric double-well potential, 

iU t + U xx + U yy - V(x, y)U - \U\ 2 U = 0, (5.7) 

where the potential V(x,y) is 

V(x, y) = -6 ( e -K*+") a -Hfl + e-K-i-B) 8 -^) , (5.8) 

which is shown in Fig. 6(a). This equation admits a family of sign-indefinite solitary waves (12.2(1 
whose power diagram is given in Fig. 6(b). It is seen that a pitchfork bifurcation occurs at the 
point 'A'. At positions 'd,e,f ' of the power diagram, profiles of the solitary waves are displayed in 
Fig. 6(d,e,f). Solitary waves at 'd,e' are anti-symmetric in x and symmetric in y, whereas the 
solitary wave at 'f is asymmetric in x and symmetric in y. Thus this pitchfork bifurcation is also 
a symmetry-breaking bifurcation. 

At this bifurcation point, we have checked numerically that zero is the largest eigenvalue of L\q, 
and (ip, LqqiP) > 0. In addition, the solitary wave at point 'A' is linearly stable, and the solitary 
waves nearby do not possess complex eigenvalues. Thus from Remark 3 and the power diagram in 
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Figure 5: Pitchfork bifurcations of solitary waves and their linear-stability behaviors in Example 3. 
(a) The power diagram (solid blue is stable and dashed red unstable); (b) profiles of the potential 
V(x) (dashed) and the solitary wave (solid) at location 'b' of the power diagram; (c) solitary 
wave profiles at locations l c\ (solid) and 'C2' (dashed) of the power diagram; (e,f,g) linear-stability 
spectra for solitary waves at locations 'b, c\, c<i respectively. 



Fig. 6(b), Theorem 4 predicts that the anti-symmetric solution branch is unstable on the left side 
of 'A' and stable on the right side of 'A', and the bifurcated asymmetric branch is stable. These 
predictions agree with our numerical stability results shown in Fig. 6(b). The numerical-stability 
results are further illustrated in Fig. 6(g,h,i), where the stability spectra for solitary waves in Fig. 
6(d,e,f) are plotted. These spectra corroborate the numerical-stability results in the power diagram 
of Fig. 6(b) and support our analytical predictions. 

In the previous examples, our comparison between analytical and numerical stability results 
was qualitative. Here for this example 4, we will also make a quantitative comparison on unstable 
eigenvalues in order to completely verify our eigenvalue formulae in Theorem 3. Specifically, we 
notice that the anti-symmetric solution branch in Fig. 6(b) is unstable on the left side of 'A', and 
this instability is induced by a positive eigenvalue which is predicated analytically by the formula 
(|3.7p in Theorem 3. Numerically we have determined this unstable eigenvalue A for various values of 
\i by the highly-accurate Newton-conjugate-gradient method [13], and these numerical eigenvalues 
are plotted in Fig. 6(c) as blue squares. Further examination of this numerical data shows that as 
fi — > fio, where fio w 1.9072149 is the propagation-constant value at the bifurcation point 'A', the 
numerical eigenvalue A behaves as 
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Figure 6: Pitchfork bifurcations of solitary waves and their linear-stability behaviors in Example 
4. (a) The two-dimensional double- well potential ()5.8j) : (b) the power diagram (solid blue is stable 
and dashed red unstable); (c) the unstable eigenvalue A versus /i on the anti-symmetric solution 
branch of (b) [squares: numerical values; solid red line: analytical prediction from formula (13. Tf) ] : 
(d,e,f) profiles of solitary waves at locations 'd,e,f of the power diagram in (b); (g,h,i) stability 
spectra for solitary waves in (d,e,f) respectively. 



where the numerical coefficient is 

cw « -0.3788744. 

In the analytical eigenvalue formula ()3.7|) . the coefficient a from formula (j3.8|) is found to be 

cw « -0.3788744. 

We see that the numerical eigenvalue formula (|5.9p and the analytical formula (|3.7p are in complete 
quantitative agreement. 

On the bifurcated asymmetric solution branch in Fig. 6(b), the solitary waves possess a pair 
of purely imaginary discrete eigenvalues which are predicted analytically by the formula (j3.9[) in 
Theorem 3. We have quantitatively compared the numerical values of those imaginary eigenvalues 
against the analytical formula (13. 9p and found complete agreement as well. Thus both eigenvalue 
formulae (|3.7p and (|3.9p in Theorem 3 are numerically verified. 
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6 Summary 



In this article, linear stability of both sign-definite (positive) and sign-indefinite solitary waves 
near pitchfork bifurcations has been analyzed for the generalized nonlinear Schrodinger equations 
(|2.ip with arbitrary forms of nonlinearity and external potentials in arbitrary spatial dimensions. 
Bifurcations of linear-stability eigenvalues associated with these pitchfork bifurcations have been 
analytically calculated, and their expressions are given by the formulae (13.71) and (13, 9ft . Based on 
these eigenvalue formulae, linear stability of solitary waves near pitchfork bifurcations is then de- 
termined. We have shown that the smooth solution branch u°(x; (j) always switches stability at the 
bifurcation point. In addition, the bifurcated solution branches tr^x;//) and the smooth branch 
have opposite (same) stability when their power slopes Pq(^o) and P±(fio) have the same (oppo- 
site) sign. One unusual feature of these pitchfork bifurcations in the GNLS equations is that the 
smooth and bifurcated solution branches can be both stable or both unstable, which contrasts such 
bifurcations in finite-dimensional dynamical systems where the smooth and bifurcated branches 
generally have opposite stability [18]. For the special case of positive solitary waves, strong and 
very explicit stability results have also been obtained. We have shown that for positive solitary 
waves, their linear stability near a pitchfork bifurcation point can be read off directly from their 
power diagram. Lastly, a number of numerical examples of pitchfork bifurcations in Eq. (|2.ip have 
been presented, and the numerical results fully support the analytical predictions both qualitatively 
and quantitatively. 
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